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Time-dependent density functional theory (TDFT) provides a way of calculating, in principle 
exactly, the linear response of interacting many-electron systems, and thus allows one to obtain 
their excitation energies. For extended systems, there exist excitations of a collective nature, such 
as bulk- and surface plasmons in metals or intersubband plasmons in doped semiconductor quantum 
wells. This paper develops a quantitatively accurate first-principles description for the frequency and 
the linewidth of such excitations in inhomogeneous weakly disordered systems. A finite linewidth 
in general has intrinsic and extrinsic sources. At low temperatures and outside the region where 
electron-phonon interaction occurs, the only intrinsic damping mechanism is provided by electron- 
electron interaction. This kind of intrinsic damping can be described within TDFT, but one needs 
to go beyond the adiabatic approximation and include retardation effects. It has been shown [G. 
Vignale, C. A. Ullrich, and S. Conti, Phys. Rev. Lett. 79, 4878 (1997)] that a density-functional 
response theory that is local in space but nonlocal in time has to be constructed in terms of the 
currents, rather than the density. This theory will be reviewed in the first part of this paper. For 
quantitatively accurate linewidths, extrinsic dissipation mechanisms, such as impurities or disorder, 
have to be included in the response theory. In the second part of this paper, we discuss how extrinsic 
dissipation can be described within the so-called memory function formalism. This formalisms will 
first be introduced and reviewed for homogeneous systems. We will then present a synthesis of 
TDFT with the memory function formalism for inhomogeneous systems, which allows one to account 
simultaneously for intrinsic and extrinsic damping of collective excitations. As an example where 
both sources of dissipation are important and where high-quality experimental data is available for 
comparison, we discuss intersubband plasmons in a 40 nm wide GaAs/Alo.3Gao.7As quantum well. 

PACS numbers: 71.15.Mb; 71.45.Gm; 73.21.Fg: 78.67.De 



I. INTRODUCTION 

The calculation of excitation energies and linewidths 
of collective excitations in extended electronic sys- 
tems is one of the outstanding problems in many- 
body thej^j^Iime-dependent density functional theory 
(TDFT)Ul3'Mffl offers a powerful and elegant approach 
to this difficult problem. To set the stage for the develop- 
ments that are to follow, we shall begin this paper with 
a summary of the key elements of TDFT. Let 



2m 



v (ri 



2 ^ 



Udn-rjl) (1) 



be the Hamiltonian of a many-electron system, where 
and pi are the canonical coordinates and momenta of 
the ith electron, m is its mass, Vo(r) is a static external 
potential, which includes contributions from randomly 
distributed impurities and other sources of disorder, and 
r/(|rj — rj|) is the Coulomb interaction potential. To cal- 



culate the excitation energies,*! one adds to H a small 
time-dependent perturbation of the form 



Hx(t) = / d 3 rui(r,£)n(r) 



(2) 



where Wi(r, t) = vi(r,uj)e~ luJt + c.c. is a periodic poten- 
tial (vi <C v ) that couples linearly to the density oper- 
ator n(r) = J2i ^( r — r i)- O ne then computes the time- 
dependent density of the system, which, in the linear 
approximation, will be given by 



n(r, t) = n {r) + m(r,w)e' 



-iut 



+ C.C. 



(3) 



where no(r) is the ground-state density, and ni(r,ui) is 
linearly related to iq(r, oj) via 



ni(r,u) = / d 3 r x(r,r',cj)vi(r',w) 



(4) 



The density-density response function x(r, r', contains 
the essential information about those excited states of 
the system that are coupled to the ground state by the 
perturbation H\{t). More specifically, in a finite system 
(atom or molecule) this response function has a discrete 
set of poles on the real frequency axis, corresponding to 
the discrete excitation energies of the system. In an ex- 
tended system, the poles merge into a continuous branch 
cut along the real axis. However, isolated poles can arise 
in the lower half of the complex frequency plane: they 
correspond to collective excitations of the system, where 
the imaginary part of the frequency defines the charac- 
teristic lifetime of the excitation. 



2 



The basic TDFT strategy for calculating x(r, t',u>) is 
to construct a noninteracting system that has the same 
ground-state density n (r), and yields the same density 
response m(r, cj) as the interacting system under study. 
The dynamics of this noninteracting system is controlled 
by an effective single-particle potential which is written 
as the sum of the total external potential vo(r) + v±(r, t) 
plus the Hartree potential 



(5) 



plus a remainder, which is known as the "exchange- 
correlation" (xc) potential v xc (r,t). It is not at all obvi- 
ous that such a potential v xcr can be constructed, but, if 
it can, then Runge and GrossEI showed that it is a unique 
functional of the time-dependent density up to within an 
additive function of time. The form of the xc potential 
depends, in general, on the initial state of the system, 
but this dependence disappears if one assumes, as we do 
here, that the system is initially in its ground state.EI The 
effective noninteracting Hamiltonian (also known as the 
Kohn-Sham Hamiltonian) that yields the exact density 
is then given by 

H KS (t) = J2 jj^ + +v H ,o(r 4 ) +f xc ,o(r*)j 

i ■* 

+ d 3 r [«i(r,t)+UH,i(r,t)+WKc,i(r,*)]ft(r) , (6) 



where both the Hartree and the xc potential have been 
written as the sum of static parts wh,o,^xc,o associ- 



ated with the ground-state density, and (small) time- 
dependent parts Uh,i,u X c,i associated with the time- 
dependent density. The static part of the Kohn-Sham 
Hamiltonian [first line of Eq. (^)] yields the exact ground- 
state density in the interacting system, while the time- 
dependent part [the second line of Eq. (^|)] yields the 
exact density response. 

In the linear response regime the xc potential can be 
written as 

v xc (r,t) = «xc,o(r)+ / dt'f d 3 r7 xc (r,r',i-t')7ii(r',t'), 



(7) 

where v xc ,o( r ) depends only on the ground-state density, 
and / xc (r,r',i — t') is the retarded xc kernel, formally 
defined as 



/ xc (r,r', *-*') = 



5v xc [n](r,t) 



Sn(r,f) 



(8) 
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Fourier transformation of both v xc and / xc with respect 
to time leads to the simpler relation 



« xc ,i(r,w) = / dV / xc (r,r', u)ni(r',u) . 



(9) 



We denote by Xks ( r , r ' , w ) the density-density response 
function of the static Kohn-Sham system [the first line 
of Eq. (g)]. The second line of the same equation can 
then be regarded as a time-dependent perturbation to the 
first. The time-dependent density response is therefore 
given by 



m(r,w)= /dVx KS (r,r>)<!ui(r>)+ I d 3 r" 



- + / xc (r',r'» 

r — r" 



n x (v",uj) 



(10) 



where the second and the third terms in the curly brackets arise from i>h,i an d w XCj i, respectively. Comparing Eqs. 
(Q) and ([To]) we obtain the following integral relation between the exact density-density response function x(r, r',cj) 
and the noninteracting response function x KS (r, r', w): 



X(r,r',u) =XKs(r,r',a;)+ AxK S (r,x,w) / d 3 y 



/ xc (x,y,w) 



(ii) 



This can also be written as 

^(r.r'.w) = X K S 1 (r,r / ,a;)-— ^--/ xc (r, r', w), (12) 

|r — r I 

where x _1 is the matrix inverse of \- 

The excitation energies are finally obtained from the 
poles of the linear response function, i.e., from the solu- 
tion of the eigenvalue problem 

|dVx- 1 (r,r',c)S(r',c)=0, (13) 



where S(r,u>) is the function that describes the spatial 
dependence of the density in the excited state. 

Eq. (|ll]) is the main formal result in the TDFT ap- 
proach to the calculation of excitation energies. From a 
fundamental point of view, calculating / xc is of course 
no easier than calculating x- The main advantage of re- 
casting linear response theory within TDFT is of a more 
practical nature: as long as the exact excitations of the 
interacting system are in qualitative correspondence to 
those of the Kohn-Sham system (a kind of "Fermi liquid" 
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assumption), the xc kernel in Eq. ( |ll| ) is expected to be 
a small correction, which can be approximated with rela- 
tive impunity.E3 The simplest approximation is to ignore 
both nonlocality in space and retardation in time. This 
leads to the widely used adiabatic local density approxi- 
mation (ALDA)til in which one poses 



de KC (n) 
dn 



(14) 



n— n(r,t) 



where e xc (n) is the xc energy densityli3 of the homoge- 
neous electron gas of density n. The right-hand side of 
Eq. (|l4j) is nothing but the local density approxima- 
tion for the ground-state xc potential evaluated at the 
time-dependent density. In terms of the xc kernel, this 
approximation implies 



/ xc (r,r',w) 



dn 2 



(15) 



which is a purely real and frequency-independent object. 

The ALDA is a remarkably successful approximation, 
despite the fact that it entirely neglects the frequency 
dependence of the xc kernel, that is, the retarded de- 
pendence of the xc potential on the density at ear- 
lier times. In atoms the ALDA has yieldedr*eason- 
ably accurate values of the excitation energies. Enl3 Most 
of the residual inaccuracy has been traced to the fact 
that the ground-state xc potential in the LDA fails at 
large distance from the— nucleus. An optimized effec- 
tive potential approach similar in spirit to Eq. (14), 



but still without retardation, yields a dramatic improve- 
ment in accuracy.Q Applications to more complex systems 
(molecules, polymers)ll3t3'l!-3 have met with similar de- 
gree of success. The essential reason seems to be that the 
frequency dependence of the xc kernel is rather weak, be- 
cause it is controlled by multi-electron excitations, which 
are either very high in energy (atoms and molecules), or 
smoothly distributed through a spectral range (extended 
systems). 

There are, however, some important features of the 
dynamical response that cannot be accounted for in any 
way by an instantaneous xc potential. Quite generally, 
the need for a dynamical theory of / xc arises in the study 
of excitations that do not have an analogue in the Kohn- 
Sham system. Perhaps the clearest example of this is 
provided by collective excitations in extended electronic 
systems, such as bulk- and surface plasmons in metals or 
inter- and intrasubband plasmons in doped semiconduc- 
tor quantum wells. In this case, the ALDA would predict 
resonance peaks of vanishing width, in glaring contradic- 
tion to experiment. 

Attempts to go beyond the ALDA to include retarda- 
tion date back to the mid-eighties. In 1985 Gross and 
Kohno proposed a dynamical local density approxima- 
tion for / xc which was designed to preserve the local re- 
lationship between w XCi i and the density, while including 
retardation in time. Their approximation reads 



satisfc-iic generalized Kohn's theorerrO ar 
rulesBE3 As a consequence, it was founded 



where f^ c (k,oj) is the xc kernel of the uniform electron 
gas calculated at the local ground-state density n (r) 
(more about / xc will be said in the next sections). Be- 
cause / xc (w) is complex, this approximation yields a finite 
linewidth for excitations that would have zero linewidth 
in the ALDA.EJl2 

Unfortunately, the Gross-Kohn approximation jl^ ) 
suffers from several inconsistencies, such, as the failure to 

211 " A . related sum 
a consequence, it was louncCi] that within 
this approximation, intersubband plasmons in quantum 
wells may become substantially overdamped. These defi- 
ciencies were ultimately traced back to the fact that a lo- 
cal approximation for the dynamical xc potential in terms 
of the density does not exist (except at to = 0, in which 
case it is the static LDA). The reason for this startling 
result is that the xc kernel of a non-homogeneous system 
is a function of infinite range in space, or, more precisely, 
the spatial Fourier transform / xc (k, k',uj) diverges when 
k — > at constant k' or vice versa. (In the homogeneous 
case one has k = k' and the-singularity disappears). 

Vignale ancLKohn (VK)c3 and Vignale, Ullrich and 
Conti (VUC)EJ showed that the nonlocality problem 
could be circumvented by working with the current den- 
sity rather than the density as a basic variable. The idea 
is to perturb the system (uh with a time-dependent vector 
potential ai(r, t) = &i(r,oj)e~ llAjt + c.c. rather than with 
a scalar potential. The perturbing Hamiltonian has the 
form 



(17) 



H x {t)= d a r ai{r,t) ■ j p (r) 



/ xc (r,r',o;) = / x l c (A ; = 0,a;)5(r-r'), 



(16) 



where j p (r) = ^ £ * [Pi<5(r - r;) + (Kr - r^p;] is the 
paramagnetic current density operator .E3 One then calcu- 
lates the current response, and determines the excitation 
energies from the poles of the current- current response 
function. 

The Kohn-Sham Hamiltonian in this time-dependent 
current density functional theory (TCDFT) contains an 
xc vector potential a XCj i(r,i), which is a (linear) func- 
tional of the full current density response ji(r, t) — 
j p i(r,t) +n (r)ai(r,i)/m: 

«xc,i,«(r,w) /xc,a/3(r,r',w)ji :/ 3(r',Lj), (18) 

J 

where / xc , 0/3(1", r', uj) is the tensorial generalization of the 
usual xc kernel (here and in the following, a, (3 denote 
Cartesian components). The static part of the Kohn- 
Sham Hamiltonian remains unchanged, and the ground- 
state density is still determined by the static xc field u XCj o- 
It turns out that the xc vector potential does admit a 
local approximation in terms of the current density: as 
we shall see in the next section, the form of this approx- 
imation is essentially determined by symmetry consid- 
erations, and can be expressed in terms of an xc stress 
tensor. □ The resulting expression for a xc 1 is local in 
space, retarded in time, satisfies the generalized Kohn's 
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and allows a consistent calculation of the 



linewidth of elementary excitations, at least the part of it 
that arises from intrinsic many-body effects. The funda- 
mental reason why all this is possible is that the relation- 
ship between the longitudinal current and the density is 
nonlocal. From the continuity equation 



dni(r,t) 
dt 



-V-ji(r,t) 



(19) 



one sees that the longitudinal component of the current 
is given by 



1 fj3 I ™l( r '^) 



(20) 



while the transverse component of the current remains 
undetermined. Thus, a local functional of ji i, will nec- 
essarily be a nonlocal functional of the density. What is 
remarkable here is that the nonlocality of func- 
tional of density can be completely eliminated by "up- 
grading" TDFT to a description in terms of a XC) i and ji. 
This is the essence of the VK and VUC theory. 

Since any scalar potential can be represented by an 
equivalent vector potential, and since the density is eas- 
ily calculated from the current [see Eq. ([l9|)l, we see 
that the ordinary TDFT is a special case of the TCDFT 
formulation.Ej An additional advantage of this formula- 
tion is that it allows one to treat the more general prob- 
lem of the response of an electronic system to an electro- 
magnetic field having both longitudinal and transverse 
components, whereas the original Runge-Gross formula- 
tion is limited to longitudinal fields, i.e., fields that can 
be expressed as the gradient of a scalar potential. 

Although the VK and VUC formulations are impor- 
tant steps enabling the calculation of the linewidth of el- 
ementary excitations in extended systems, they are still 
not sufficient to achieve quantitative accuracy in cases 
of practical interest. For example, the calculation of the 
linewidth of the intersubband (ISB) plasmon in a 40- 
nm GaAs/Al .3Gao.7As quantum well reported in Ref. 
p3| , based on VUC formalism, yielded a linewidth about 
5 times smaller than the experimental value. The rea- 
son for this disappointing result is that the theory, as 
it stands, does not take into account other intrinsic and 
extrinsic sources of damping, such as electron-phonon in- 
teractions, electron-impurity scattering, and, in the case 
of quantum wells, interfacial roughness. All these inter- 
actions contribute to the linewidth and must be included 
in any calculation that aspires to quantitative accuracy. 

In this paper we take a first step in this direction by 
showing how two of the most prominent contributions 
to the low-temperature linewidth of plasmons in quan- 
tum wells, namely electron-impurity scattering and in- 
terfacial roughness, can be built into the current-density 
functional formalism. 

Our approach is based on a the "marriage" of the 
TCDFT formalism with the memory function formal- 
ism described, for example, by Forster.c3 In the homo- 
geneous electron gas limit this approach reduces to the 



Belitz-Das Sarmac3 treatment of the effect of impuri- 
ties on bulk plasmons, which, in turn, can be viewed as 
the high-frequency extension of the Mermin relaxation- 
time approximationca for the density-density response 
function of an electron gas in the presence of ran- 
domly distributed impurities. Of course, our interest 
lies in strongly inhomogeneous systems, such as quan- 
tum wells,E2IEil which exhibit the intersubband plasmon 
resonance. Such resonances are of practical interest in 
connection with the design of infrared detector devices. 

Our strategy is to derive an integral equation which 
relates the current response function of the disordered 
interacting many-electron system to that of the same sys- 
tem in the absence of disorder: the latter is calculated 
by the standard TCDFT outlined above. We shall show 
that this approach (despite some inevitable approxima- 
tions in the treatment of disorder) meets with consider- 
able success: the linewidth of the intersubband plasmon 
is considerably enhanced by disorder - in particular, by 
interfacial roughness - and agrees quantitatively with the 
measured one. More importantly, the qualitative behav- 
ior of the linewidth as a function of an external electric 
field that controls the shape of the quantum well is cor- 
rectly reproduced. 

The remainder of this paper is organized as follows: In 
Section [0] we review the main aspects of the TCDFT for- 
malism for the linear response of many-electron systems. 
In Section [iFJ we review the memory function formalism 
and demonstrate its application to the case of a homo- 
geneous electron gas with randomly distributed impuri- 
ties. These sections are meant to make the paper self- 
contained and to provide the necessary background for 
the following more technical parts. In Section IV we com- 
bine the memory function formalism with TCDFT for 
inhomogeneous systems and derive the key integral equa- 
tions for the current-current response functions. Finally, 
m Section [V] we demonstrate the power of the method by 
calculating the linewidth of the intersubband plasmon in 
a quantum well and comparing to recent experimental 
results, and in Section [Vl] we give our conclusions. 



II. TCDFT BEYOND THE ADIABATIC LDA 

A. Exchange-correlation kernels in the 
homogeneous electron gas 

For orientation, let us first consider the xc kernels of a 
homogeneous electron gas. Because of translational in- 
variance, it is convenient to work with the Fourier trans- 
form j(k, oj) of the current density. The linear response of 
this quantity to a vector potential ai (k, to) can be written 
as 



ji,a(k,w) 



^Xa/3(k,w)ai !/3 (k,a;) 



(21) 



where Xa/3(k, w) is the current-current response tensor. 
Due to rotational invariance, the responses of the longi- 
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tudinal (parallel to k) and transverse (perpendicular to 
k) components of the current are completely indepen- 
dent, and one can write 



Ji,i(T)(k,o;) = XL(T)(k, w)ai iL(T )(k, w) , 



(22) 



where L(T) denotes the longitudinal (transverse) com- 
ponent, and XL(T)(k, ui) is the longitudinal (transverse) 
response function.—. According to the general linear re- 
sponse formalism,L3 Xl(t) is given by 



Xi(T)(k,w) 



J2\(l\kL(T)M\0)\ 2 



1 



1 



U! - LUX) + IT] Ul + LU10 + IT] 



,(23) 



where |0) is the ground state, |Z) is the l-th excited state, 
and loio is the excitation energy Ei — Eq. The Fourier 
transform of the paramagnetic current operator j p (k) is 
given by 



jp( k ) = 77- V] [Pi 



(24) 



Note that the total current response is the sum of the 
"London current" nai/m and the paramagnetic current 
[the expectation value of (Q)]. A key feature of Eq. (|23] ) 
is that in the limit of k —* and finite uj it has the limiting 
form 



- a L (u>)k 2 



XT(k,ui) — » harHfc 

m 



(25) 



where aifr)(w) are functions of frequency only, e(u>) — 
1 — limfc_,.o nv(k)k 2 /mui 2 is the homogeneous dielec- 
tric function, and v(k) is the Fourier transform of the 
Coulomb interaction. Notice that the difference between 
the longitudinal and transverse results at k = is due to 
the long range of the Coulomb interaction: the difference 
vanishes if v(k) diverges more slowly than 1/fc 2 . 

Translational invariance is the essential reason for the 
small-fc behavior of the %'s. In the k — ► limit the cur- 
rent operator reduces to the total momentum operator 
plus a correction that vanishes linearly with k. Thus the 
first term on the right-hand side of Eq. ( poj ) is the re- 
sponse of the center-of-mass momentum, which obeys a 
simple equation of motion under the action of the exter- 
nal force, while the second term, of order k 2 , comes from 
the residual part of the operator, which is linear in k. 
There are no cross terms, since the dynamics of the cen- 
ter of mass is decoupled from that of the internal degrees 
of freedom. 

Let us now turn to the xc potentials. The idea is to 
express the exact current response ( p2| ) as the response 
of a noninteracting electron gas to an effective vector 
potential, written as 

aeff,i(k,w) = a x (k,w) + a H ,i(k,w) + a XCj i(k,w) . (26) 



The Hartree component ap^i is purely longitudinal (since 
it is just another way of describing the scalar Hartree 
potential) and is given by 



a H ,i(k,w) = — v(k)ji jL (k,u)k 



(27) 



The xc potential can be decomposed into its longitudinal 
and transverse components (with respect to the direction 
of k) as follows: 



a xc ,i(k,w) = — /4 L (fc,w)ji,i(k,a;)k 
ui z L 

+ /^ r (*,w)ji,r(k,«)]. (28) 



The factor k 2 /uj 2 has been introduced, in analogy to (|27 
so that the longitudinal component of the xc vector po- 
tential is equivalent to the scalar xc potential 



L (k,w) 



WO xc ,l,L 



f* L (k,w) ni Qt,w) (29) 



(ii = niui/k), thus making f£ c L (k,uj) identical with the 
usual /x C (&! uj) of the ordinary density functional theory.i 
In analogy to Eq. (p"2]), the relationship between the in- 
teracting current-current response function and its non- 
interacting counterpart % K s takes the form 



X L (V)( fc ' 



^ks,L(T) 
k 2 



(k,uj) 

VL(T)(k) + fx cMT) (k,uj) 



(30) 



where we have defined Vj J (k) = v(k) and fy(fe) = 0. 

Notice that, according to Eqs. (p5|), x^ 1 (fc,cj) and 
XZ S L (k, uj) — k 2 VL{k)/ui 2 have the same limit me(uj)/n + 
0(k 2 ) for k — > and finite uj. Similarly, (A;, w) ano - 
X~g T (k,u)) — k 2 VT(k)/u> 2 have the same limit m/n + 
0{k 2 ). Thus, we see that Eq. ( j30| ) is consistent with 
the limiting forms ( p5[ ) if and only if the k — > limits of 
the xc kernels f£ c L ^ are finite functions of frequency: 



fe^O 



^/x\l ( t)M- (31) 



Because of the central role these functions play in the 
developments to follow, we now describe their properties 
in detail. 



B. Properties of the homogeneous xc kernels 

The calculation of the xc kernels /xcL(T)( tJ ) ^ s a 
very difficult problem in many-body theory. Approx- 
imate calculations have been done using (i) Interpo- 
lation-, schemes between exact, high and low frequency 
limitsa (ii) Perturbation theoryEj-EJ (iii) Mode decoupling 
approximations.^ Here we simply summarize the main 
results that have been established to date, and refer the 
reader to the original references. 
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1. The high-frequency limit is a purely real constant 
given by 

lim /xc,L(T)M = ^[ d MT)((ke) - (fce)o) + e L{T) (pe)] , 

which is also known as the third moment sum rule, (fee) 
and (pe) are the expectation values of the kinetic and 
potential energy, respectively, and (fce)o is the nonin- 
teracting kinetic energy. In three dimensions, d^ = 4, 
e l = 8/15, dr = 4/3, er = —4/15. In two dimensions, 



d L = 6, e L = 5/4, d T = 2,e T = -1/4 (see Ref. |35|). 

The behavior of the imaginary part of the longitudi-| 
nal xc kernel was first determined by Glick and Longed 
(3D) and Holas and SingwEJ (2D), making use of second- 
order perturbation theory, which becomes, in all likeli- 
hood, exact in the high-frequency limit. More recently, 
their calculation has been confirmed and extended to the 
transverse kerne£3 by a different method based on the 
equations of motion for the current response function. 
The result is 



lim Im/. 



xc,L(T) 



4-D--D/2 
-CLL.T 7T LO ' - 



X-D 



(33) 



where a = Ti 2 /me 2 is the Bohr radius, lo — fiLoao/e 2 
is the dimensionless frequency, and D is the number of 
spatial dimensions. The coefficients are az, = 23/30 and 
a,T = 16/15 in three dimensions, and = 11/32 and 
ay = 9/32 in two dimensions. 

2. The-tj — > limit was first worked out by Conti and 
Vignale,t3 and is subtly different from the static limit, 
which is obtained by setting lo = before letting k — > 0. 
The result for / xc £ in D dimensions is 



lim Jxc,i(w) = 4'cW + 



2(D 



D 



■^lim/jLrH, (34) 



where e xc (n) = d 2 e xc (n) / dn 2 . The second term on the 
right-hand side of this expression is proportional to 



in three dimensions, and 



2E F {3F 2 - $F X 



25?i 



Fx 



Ep f F2 — F% 



2n 



Fx 



(35) 



(36) 



in two dimensions, where Fi are the-conventional Landau 
parameters of the electron liquid.E3 Note that the first 
term on the right-hand side of Eq. (|34|) is the usual com- 
pressibility obtained from the static limit lim^o lim^^o 
/^ cL (fc,o;). Thus, Eq. (34) vividly shows the non- 
commutativity of the k — > and lo — > limits. 

3. Due to the causality properties of the linear response 
functions, the xc kernels must be analytic functions of lo 
in the upper half of the complex plane. This leads to the 
Kramers-Kronig relations, which relate the real parts of 
the xc kernels to their imaginary parts: 



Re/* (w) = /*(c 



7T 



LO 



(37) 



where P denotes the "principal part!! integral. 

4. In a very recent developmental! the low- frequency 
behavior of the imaginary parts of the xc kernels has also 
been calculated exactly to leading order in the strength 
of the Coulomb interaction. The results are 



Im/xc,r(w) 



1 



n 



( na 



D\2 



(38) 



and 



Im/ XC , L M -> 2( ^ n ^ Tmf xc , T (w), (39) 



D 



where the dimensionless "shear viscosity' 
by 



55 Tj_ 



jyr is given 



(Voo) 



45tt 3 \ 
2 . 



A + j- ) tan 1 A 



— sin ; 

a VTTy 



A\/2TA2 



tan 



\V2T~X 2 



40) 



in three dimensions, and 



12tt 2 



ln(A + 1) - 



1 + A 



dx- 



\ 2 x 



(Xx + l)(\VT^ 2 + 1) 



(41) 



in two dimensions. In the above expressions A is defined 
as A = 2kp/k s , where k s is the screening wavevector, e.g., 
in RPA, k s — (ATrkpao) 1 ^ 2 /t^olq in three dimensions, and 
k s = 2/ciq in two dimensions. The derivation of these 
results is presented in Ref. [}7[ 

5. Parametrized expressions. To keep our presentation 
self-contained, we also include the explicit parametriza- 
tion for / xc ,l(w) that has been used in the calculations 
of Section V. Thja-,is the original Iwamoto-Gross-Kohn 
parametrizationJ3H and has the form 



Im/ Xc ,z,(w) = 



a{n)LO 



[l + 6(n)w 2 ] 5 /4 



(42) 



with the coefficients a(n) and b(n) determined by the 
compressibility and third moment sum rules, and the 
Kramers-Kronig dispersion relations. The real part of 
/xc,l(^) is then calculated with the help of the disper- 
sion relation (|3^). More recent analytic expressions for 
/xc,l (w) .and / xc ,t(w) have been obtained-by Nifosi, Conti 
and TosiEa and Qian and Vignale (QV).E3 These new ex- 
pressions possess considerable structure in the frequency 
dependence due to two-plasmon excitations. The QV ex- 
pression reproduces the exact perturbative limit of / xc (w) 
in the limit lo — > 0. Additional details about these ex- 
pressions can be found in the original references. 



7 



C. The exchange-correlation field for a homogeneous electron gas 



As a preparation for the study of inhomogeneous systems let us now examine the real-space form of the xc vector 
potential. It is convenient for this purpose to introduce the xc electric field 

E xc ,i(k,o;) = iwa xc ,i(k,w) 

= -^{k[k-j 1 (k I a;)]/* )1 (w) + fc 2 /^,r(w)j 1 ,T(k ) w)}. (43) 

Splitting off the familiar ALDA contribution 

E4 L 1 DA (k,c)= l k4' c ^, (44) 

UJ 

we can write 

E XCil (k,c) = E£f A (k, W ) 

- ^{klk.j^k^I/^H-e^+fcVlTHji^lkc,)}. (45) 

Introducing at this point the velocity field 

n=' i - (46) 
n 

and Fourier transforming Eq. ( [45| ) to real space, we get 
E xc>1 (r,o,) = E^? A (r, W ) 

+ £{ [/*. £ (o;) - / X \ X M - 4' c ] V(V • u) + /4 iT (c)V 2 u}. (47) 

It is easy to verify that this expression can be re-written as 

^xc : i,a(r, uj) - i/ XCj i, Q (r, uj) + - ' ( 48 J 

where the xc stress tensor a XCia /3 is defined as 

-, s ( du a dun 2^ \ ~, 

o- X c, a/ 9 = I + -r^- - —V • ud a pj + C(n,w)V • ud Q/3 , (49) 



and 



C(n,w) 



fj{n,uj) 



— / X c,tH 



(50) 



2(D-1) 



iu) \ : ' ' D 

(51) 

are generalized (i.e., frequency-dependent and complex) 
visco-elastic constants of the electron liquid.E3 In partic- 
ular, the real parts of f] and Q [related to the imaginary 
parts of f£ c L ^ T j (uj) / uj] play the role of shear and bulk vis- 
cosities, respectively, while the imaginary parts of wry and 
uj( [related to the real parts of f£ c L , T ^ (uj)] are interpreted 
as post- ALDA xc contributions to frequency-dependent 
elastic constants /i (shear modulus) and K (bulk mod- 
ulus) of the electron liquid: /idyn = "^/xcT^) an< ^ 



K dyn = n 2 Re{/ x \ T M - [2(D-l)/D]f^ c T (uj) - (see 
Ref. Bq for details; the full elastic constants fi and K 
given there also include the kinetic and the ALDA part 
of the xc contribution). 

Because Im/ xc .l(t)(^) vanish linearly for uj — > (point 
4 of Section II. B), we see that the viscosity coefficients 
fj and C stay finite in the uj — > limit. Eqs. 
of Section II. B imply that the ui — ► limit of the 
viscosity lim^^o C( w ) is exactly zero at least to within the 
accuracy of our pcrturbative calculation. By virtue of the 
limiting form (|34|), we also see that lim u _»o K& yn (uj) = 0, 
implying that the bulk modulus of the electron liquid is 
entirely accounted for by the ALDA contribution e xc . 

The fact that the shear modulus does not vanish for 
ui — ► but tends to the finite value of Eq. ( |34| ) is perhaps 
surprising. One ordinarily thinks of liquids as having zero 



|3£ 
bulk 
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shear modulus. The reason for this strange behavior is 
that we are taking the k — > limit before the uj — > 
limit. Thus, the system remains "dynamical" down to 
zero frequency. Of course, this would not be true if the 
uj — > limit were taken at hnite k. In that limit, f£ T is 
no longer related to the shear modulus, but to the static 
diamagnetic susceptibility, which is extremely small. The 
truly static shear modulus is zero, as expected. 



D. The exchange-correlation held in the 
inhomogeneous electron gas 

The main result of the previous section, Eq. is 
artfully written so that it can immediately be turned into 
a local density approximation for the xc electric field of an 
inhomogeneous electron liquid through the replacement 
ri — > n (r), where n (r) is the ground-state density of 
the inhomogeneous liquid. Of course, the xc kernels must 
also be evaluated at the local density. 

An important question is this: Why should the re- 
placement n — > no (r) be done in Eq. ( Eq ) rather than in 
one of the many equivalent expressions one can generate 
starting from Eq. (plEf)? For example, why not write the 
second term on the right-hand side of Eq. (|47|) in the 
equivalent form 

^{ [/x*UM - /ItM - 4] V(V-j 1 ) + /l T (.;)V 2 j 1 } 

(52) 

before substituting n by «o( r )? The answer is that this 
and similar ambiguities are completely removed by gen- 
eral physical requirements which we now discuss. 

First of all, because the Coulomb interaction obeys 
Newton's third law, the net force exerted by the xc elec- 
tric field on the system must vanish. At the local level, 
Newton's third law implies that a small volume of the 
electron liquid cannot exert a net force on itself. Accord- 
ingly, the net force acting on an arbitrary volume ele- 
ment must be expressible as the integral of the external 



stresses exerted by the surrounding fluid on the surface 
of the volume element. The mathematical expression of 
this requirement is that the force density must be the 
divergence of a local stress tensor as in Eq. (Q). 

A similar argument can be applied to the net torque 
acting on a volume element of the fluid. Again, this must 
be expressible in terms of a surface integral, and it is not 
difficult to see that the condition for this to happen, is 
that the stress tensor be a symmetric rank-2 ten.sor.E3 

Finally, Galilean invariance requires the stress tensor 
to vanish identically when the fluid moves as a whole, i.e., 
when the velocity field is spatially uniform. It is for this 
reason that the stress tensor must contain derivatives of 
the velocity field and not of the current. 

These criteria unambiguously establish Eq. ([f8j) as the 
correct expression in which the substitution n — > n>o(r) 
should be made. This expression was originally derived 
by Vignale and Kohn (VK) by a more laborious and ap- 
parently quite different path, which more clearly exposed 
the underlying approximations and the conditions for 
their validity. VK considered a weakly inhomogeneous 
electron liquid modulated by a charge-density wave of 
small amplitude 7 and small wave vector q. Both k and 
q were assumed to be small not only relative to the Fermi 
momentum kp but also relative to ujjup (uf being the 
Fermi velocity). The latter condition assures that the 
phase velocity of the density disturbance is much faster 
than the Fermi velocity, so that no form of static screen- 
ing can occur. Under these assumptions, all the compo- 
nents of the tensorial kernel / xc , Q /3 could be calculated, 
up to first order in the amplitude of the charge density 
wave, and to second order in the wave vectors k and q. 
The calculations were greatly facilitated by a set of sum 
rules which are mathematically equivalent to the zero- 
force and zero-torque requirements discussed above. The 
result of the analysis was that the diagonal matrix ele- 
ments / xc ,Q/3(k, k, uj) remain equal to f^ c {k, uj) to first or- 
der in 7, but the off-diagonal elements / xc ,a/3(k + q, k, uj) 
acquire a finite value, given byE3 



/xc,a/?(k + q,k,w) 



f) fh 

(SfZcL - fL T )q a q f 3 + f^ T q 2 s a p - n -gSZl k • (k + d)5 af} 



J_ 

UJ 2 

+ A(n, uj)(k a + q a )kp - B[ u.x)k,A k < + .7 ,) \. 

r 



(53) 



where 6f£ cL = f£ cL (uJ,n) - e'Un), A(n,w) = 
[n(2d% cT /dn - dftJdn + 3f* cT - 5f* cL ] and B(n,uj) = 
[ndf£ cT /dn+3f£ cT — 5f£ cL \. A remarkable feature of this 
result is that the off-diagonal matrix elements of / xc ,a/3 
do not exhibit any singularity for k or q tending to zero 
in any order. This is in marked contrast with the off- 
diagonal elements of the scalar (density) xc kernel which, 



when calculated for the same system, exhibit a power sin- 
gularity of the form k ■ q/fc 2 for k — > at finite q. This is 
the fundamental reason why the dynamical local density 
approximation is possible in terms of the current, but not 
in terms of the density. 

Eq. (|53|) can be translated into a real space expression 
for E XCj i(r,cj). More details of the derivation, which is 



9 



quite laborious, are given in Ref J22|. Finally, the result- 
ing expression can be rearrangedcJm the elegant form of 
Eq. fl48|) . The conditions of validity of the real-space ap- 
proach are |Vno(r)|/no(r) much smaller than fcp(r) and 
(jj/uf(j), where fcp(r) and uf(t) are the local Fermi mo- 
mentum and velocity. In Ref. E3l the practical relevance 
of these conditions was investigated in detail. It was 
found that the approach could be successfully applied to 
describe intersubband plasmons in wide single quantum 
wells, but failed for narrow double quantum wells. The 
failure in the latter case was traced back to a strong vi- 
olation of the above criteria of validity in the region of 
the barrier between the two wells. The physical reason 
is that electronic motion through tunnelling barriers im- 
plies strong internal compression of the electron liquid, 
which locally destroys coherence of the electron dynamics 
and leads to a breakdown of the simple hydrodynamical 
picture. In such a situation, a hybrid between the VUC 
and the more robust Gross-Kohn approximation (13) pro- 
vides a pragmatic and practically useful remedyEj 

Because the occurrence of two spatial derivatives of 
the velocity field in the post-ALDA term is dictated by 
general principles, Eq. ( 4§|) is expected to remain valid 
even for large values of u, i.e., in the nonlinear regime, 
provided that u and n are sufficiently slowly varying. The 
argument goes as follows. Suppose we tried to extend 
Eq. ( |48| ) into the nonlinear regime by including terms 
of order u 2 . Because the stress tensor must depend on 
first derivatives of it, such corrections would have to go as 
(Vm) 2 . But then the force density, given by the derivative 
of the stress tensor, would have to involve at least three 
derivatives. Thus, for sufficiently small spatial variation 
of the density and velocity fields, the nonlinear terms can 
be neglected. 

Since the ALDA is an intrinsically nonlinear approx- 
imation, VUC proposed that Eq. (f|8|), written in the 
time domain, could provide an appropriate description 
of both linear and nonlinear response properties. A non- 
linear, retarded expression for f xc ,i was also proposed 
by Dobson et aZ.c2l The two approximations coincide in 
"one-dimensional systems" (i.e., when one has a unidi- 
rectional current density field that depends only on one 
coordinate), but differ in the general case. 



III. MEMORY FUNCTION FORMALISM AND 
TCDFT FOR HOMOGENEOUS SYSTEMS 

In the preceding sections, we outlined a linear response 
formalism within TCDFT that goes beyond the adia- 
batic approximation and allows one to account for in- 
trinsic damping of collective excitations in electronic sys- 
tems, caused by dynamical many-body effects. As men- 
tioned in the Introduction, this is usually not sufficient 
to achieve quantitative agreement with experimentally 
measured linewidths. In reality, intrinsic damping is of- 
ten overshadowed by strong extrinsic dissipation mecha- 
nisms, such as impurities or disorder (in this paper, we 



consider the low-temperature case only and limit the dis- 
cussion to systems where LO phonon scattering does not 
occur). 

Effects of impurities and disorder in the linear dy- 
namics of a many-electron system are conveniently dis- 
cussed in the language of relaxation functions J-H which 
then naturally, leads to the so-called memory function 
formalism YWA In this paper, we perform a conceptually 
new step and unite the memory function formalism with 
TCDFT in the linear response regime, which will then al- 
low us to treat both intrinsic and extrinsic damping from 
first principles and on an equal footing. This is necessary 
for an accurate description of experiments performed on 
very clean samples (such as the quantum well we shall dis- 
cuss in section ^) , where intrinsic and extrinsic damping 
may be of comparable magnitude. 

The purpose of this section is twofold: To make this 
paper self-contained, we first review the memory func- 
tion formalism for the homogeneous case. We then 
make contact between relaxation and response functions, 
thereby integrating the memory function formalism with 
TCDFT. In Section IV we shall extend the approach to 
systems that are inhomogeneous in one spatial direction 
(such as quantum wells), and show how it can be applied 
to discuss intrinsic and extrinsic damping of collective 
charge-density excitations in such systems. 



A. Relaxation and Linear Response 

Suppose we are interested in the dynamics of a set of 
N + 1 observables, {A (r,t), . . . , A^ir^t)}, each Ai(r, t) 
coupling to a small perturbing external field Scii(r,t). 
The Hamiltonian in the presence of these fields is given 
(in Schrodinger representation) by 



H(t) =Ho + ^2 J d3r Mr)Sa t {r, t) 



(54) 



Consider now the case where the external fields are adi- 
abatically turned on beginning at t = — oo, and then 
abruptly switched off at t = 0: 



Sdi(r, t) 



8<n{v) e"* 




till (») 



The system starts out from non-equilibrium at t = + 
and, being left to itself, relaxes back towards equilibrium. 
The first-order change of the non-equilibrium osnefjiatipn 
value of an observable can then be written asEZrclrEdej 



S(Mr,t)) 



noncq 



= E J^r'C ij {T,T',t)5a j {v') , (56) 



where the correlation (or Kubo) function in the presence 
of disorder is defined as 



Cy(r,r',f) 



d/3' (Mr,t)A 3 (r',-i(3')) c 
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(A i (r,t)) eq .(^(r',-0 , ))eq.J j (57) 

/ disorder 



This may be rewritten as 

C ii (r,r',t) = (i i (r)|e-^ t |i i (r') 



(58) 



where £ is the Liouvillc operator governing time evolu- 
tion of the system via 



Ai(t)=iLAi{t) = [A^.HyiiK), 



(59) 



and the scalar product (. . . | . . .) is defined by Eq. (57). 
We impose normalization on the set of variables, i.e. 
(Ai(r)\Aj(v')) — Sij. In the following, we are interested 
in the zero temperature limit (8 — ► oo). 

The correlation function IpTtj is related to the absorp- 
tive part of the usual quantum mechanical response func- 
tions Xij( r , r \t) as follows: 



id t C i3 (v,v',t)=^!Av,v',t). 



(60) 



Since relaxation occurs for t > 0, it is customary to intro- 
duce the Laplace transform of the quantities of interest, 
which, e.g., for the correlation function is defined as 



/>oo 

Cy(r,r',£)= / dfce«'Ci i (r,r , ,i), 
Jo 



(61) 



where £ is a complex number in the upper half of the 
complex plane. One then finds the following relationship 
between the Kubo relaxation functions and the response 
functions: 

^•(ry.f) = { Xij (T,T',0-Xij^^'^0))/(iO . (62) 

For the remainder of this section, we assume that the 
system is spatially homogeneous (a more general case will 
be considered in Sec. IV). Eq. ( pq ) can then be Fourier 
transformed into momentum space, and one obtains 



noncq 



= ^Cy(q,0 S aj {q) 



(63) 



The Laplace transform of the correlation function (f58|) is 
then given by 



C tJ (q,0 = (i 4 (q 



AM 



(64) 



B. Projectors and Memory Functions 

The observables {A (q), . . . , Ajy(q)} can be regarded 
as vectors in a Hilbert space. The Liouvillian C acts 
as a linear operator in that space, see Eq. (|59|). We 
define a projection operator V onto the space spanned 
by {A)(q), . . .,A N (q)} as 



P = ^|i l ( q ))(i l (q)|, 



and its complement Q = 1 — V projects perpendicular to 
it. One can then formally write Eq. ( |64| ) as 



t-CQ-CV 



AM 



(66) 



Following Forster,0 o ne p erforms a few straightforward 
manipulations in Eq. ( p6[ ) and finds 

CiMO = | % + 1 E {<ii(q)|£|ifc(q)> 



Q 



i 



QCQ-t 



Q 



AM))>c k M^)-m 



Defining 



C«(q,0 = iCiM,0 
nifc(q) = (i i (q)|£|i fc (q)) 

M lfc (q,£) = (AM ' 



Q 



(68) 
(69) 

-l,.(q)) (70) 



we rewrite Eq. (^) as follows: 

E {to* - M<0 + MiM,0}c k j(<i,0 = -Sij . (7i) 

k 

f2ij(q) can be viewed as characteristic frequency or 
restoring force matrix of the system. Our particular 
interest, however, lies in the memory function matrix 
Mifc(q, £), which introduces dissipation into the electron 
dynamics of the system. In general, dissipation origi- 
nates from intrinsic as well as extrinsic scattering mech- 
anisms. The former, caused by electron-electron inter- 
action alone, are present even in a perfectly "clean" 
system.cj In the previous section we discussed the treat- 
ment of intrinsic dissipation in the framework of TCDFT. 
Now we describe how additional extrinsic dissipation (e.g. 
caused by scattering off disorder or charged impurities) 
can be taken into account simultaneously. 

Mik has the formal structure of a correlation func- 
tion between two projected forces, QAi{q) and QA,(q). 
These forces act perpendicular to the vector space of 
variables {^4i}, thus providing a coupling to other de- 
grees of freedom of the system (which effectively form a 
"thermal bath" ) . Accordingly, the frequency dynamics 
of Mjfc(q, £) is determined by Q£Q, where those fluctu- 
ations of the Liouville operator are projected out that 
occur only within the space of variables {^4i}, therefore 
describing the internal dynamics of the "bath" . 

The correlation functions CV,(q, £) arc determined by 
a set of (N + l) 2 coupled equations, Eq. (71), whose 



solution will be discussed for an example in Sec. Ill C . 
The observables Ai are not restricted to be scalars, but 
can also be vectors (or nth-rank tensors). In general, all 
correlation functions CV,'(A:, to) 



as well as f2jj and My 



are 



(g5) tensors whose rank equals the sum of the ranks of Ai and 
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Assuming that explicit solutions for the Cij have been 
found, the final step is then to make contact with lin- 
ear response theory, which involves Fourier transforms 
(with frequency to) rather than Laplace transforms (with 
frequency £) of the associated response and correlation 
functions. Fortunately, the relationship between Fourier 
and Laplace transforms is a straightforward linear one, 
so that the response functions Xij(q, w ) can be simply 
obtained from 



Cij (q, w) = [xij (q, cj) - Xij (q, 0)] /u 



(72) 



Here, u> are real frequencies, and Cy(q, uS) = Cy(q, £ 
uj + iO+). 



Eq. (|7q) , with i,j = 0,1, represents a system of 4 
equations coupling the 4 possible longitudinal correlation 
functions (p-p, p-j L , j L -p and j L -j L )- Since there is no 
coupling between density and transverse currents, there 



is only a single transverse correlation function, j -j 



T -T 



determined by Eq. 



Using the continuity equation, 



43(q) = -?i L (q) , 



(77) 



we convince ourselves that indeed (Ao|A^) = 0. Further- 
more, since CAq is proportional to A\ , the first compo- 
nent of the "perpendicular" force QAo is identically zero, 
so that 



C. Generalized Relaxation Time Approximation 

In the following, we discuss the case where there are 
only two observables of interest: density fluctuations /5(q) 
and current density j(q). Including normalization, we 
have 



^o(q) = 



p(q) 



</S(q)l/Kq)> 1/2 



and 



^i(q) = 



j(q) 



(73) 



(74) 



where n is the uniform density of the system, and 
(/3(q)|/5(q)) = %(q, 0), the static density-density response 
function (in the presence of disorder). Eq. ( |7l| ) then de- 
scribes the four associated correlation functions Coo, Cm, 
Cio and Cii, i.e. p-p, p-j, j-p and j-j. 

Since we deal with the homogeneous and isotropic case, 
Eq. (fn]) decouples into two independent equations for 
the longitudinal (L) and transverse (T) components of 
the correlation functions: 



l 

E 

k=0 



q) + M^(q,o;; 



fi£ (q) + Ml (q, to) C£ (q, w) = - 1 



(75) 



(76) 



±\±r\ 



xro(q > u)=M o i 1 (<fcw)=M&(q,w)=0. (78) 
In the following, we will be concerned with the limit of 
weak disorder. In this limit, it is a good approximation to 
assume that all static correlation functions are unaffected 
by disorder. This means that £l L ( T \q) contains effects 
of Coulomb interaction only. Likewise, we assume 



X 



L(T) 



(q,o) 



*r (Tc) < 



q,0) 



(79) 



where the superscript "c" denotes the "clean" response 
function. In general, static disorder effects (mainly 
changes of the density of state at the Fermi surface) are 
weak, provided (tef)^ 1 <C 1, where £f is the Fermi en- 
ergy, and r is a characteristic disorder scattering time for 
the system under study. 

In general, the memory functions contain both intrinsic 
(Coulomb interaction) and extrinsic (disorder) contribu- 
tions. In the limit of weak disorder, one can separate 
them as follows: 



M 



L(T) 
11 



(q,w) 



M in y ; (q,w) 



■MWfaw). (80) 



We combine the intrinsic part with fljfc, thereby defining 
a frequency-dependent, dissipative restoring force matrix 
(q, u>), which contains effects of Coulomb inter- 



n 



ik 



MT) 



actions only: 

^; L(T) (q,c) = ttf fe (T) (q) + A^ (T) (q, u)8 ix 5 kl . (81) 
One then obtains from Eqs. (|75|) and (^6h: 



i 

J2 {[w + M c L x (q, u)] S lk - ^ L (q, w)}tf&.(q, u) = -S tJ + 8 iQ M^ u)C± (q, w) , i,j = 0, 1 



k=Q 



(82) 



lo + Af£(q, u) - tti{ T (q, wj\ Cf x (q, cj) = -1 . (83) and the subscript"ex" of and in the following. 

To solve Eqs. (g2J) and (|8^), we introduce "clean" lon- 
For notational brevity, we suppress the (q, w) dependence gitudinal (transverse) reference functions C^ c(Tc) (q, w), 
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defined in the absence of dissipation (i.e., M L = 
M T ), as follows: 



2[w*ifc-n« ,i (q,w)]^/(q 1 w) 



fc=0 



hi , i,3 = 0, 1 , 
(84) 



w - ^^(q.w) C{f(q,w) = -1 



(85) 



The desired correlation functions are then expressed in 
terms of these reference functions. In the longitudinal 
case, we find 



C£(q,w) 
or explicitly 

^(q,w) 
C&(q,w) 



C£ c (q, w + M L ) - Cf c (q, w + M L ) 



xM 1 C£.(q,«), i,j = 0,l, 



°00 



1 + M^Cgf 



CLc 
10 



1 + M^Cft 

r~iLc 

1 + M^Cgf 



U n j- Oil ^ ^QQ ~ Wo «^( 

l + M^ c 



(86) 

(87) 
(88) 
(89) 
01 ,(90) 



where all "clean" functions carry the arguments (q, u> + 
M L ). For the transverse case, on the other hand, we 
simply obtain 



C£(q,w)=C£ c (q,W + M T ) 



(91) 



We now make contact with the response functions as 
follows. For the case of density-density response, we have 
fromEq. (fr§): 

C L (q, w) = - (q, w) - Xoo 0)1 / x£> (q, 0) , (92) 
where the function Xoo(q'O) m the denominator arises 



from the normalization of the variable Aq, see Eq. (73). 
Similar expressions are obtained for the other longitudi- 
nal and transverse correlation functions. We then find 
fromEq. @: 



1 



1 



xhM^) " (w + m l ) x^(q^ + m l ) 

M L 1 



(93) 



Eq. (93) is formally in agreement with Belitz and Das 
Sarma[equation (2.3) in Ref. j2j|. However, Xoo(q||, a; ) 
here denotes the exact, fully interacting response func- 
tion of the homogeneous electron gas, not just the RPA 
response function. 



^From (|93|), one easily derives explicit expressions 
for longitudinal density-current and current-current re- 
sponse functions using 



xob(q^) = -xfo(q> w ) = -xoi(q, w ) = -^xn(q^)- 

Lu UJ U) 

(94) 

It is not difficult to show that the results for Xioi Xoi 
and X r\ o btained in this fashion are consistent with Eqs. 
©-©■ 

In the same way one finds the transverse current- 
current response function from Eq. (|9l|): 



I 2 ,.L 



xii(q> w ) 



(w + M T ) 
M T 



M T ) 



Xi ?(q,w + M T ) 
X^(q,0). (95) 



Expressions that are formally similar to Eqs. i-(P3j)-(|95|) 
were recently derived by Conti and VignaleESin thft 
framework of Mermin's relaxation time approximation.E3 
In this formalism, the role of M(q,oj) is taken by a 
frequency- and momentum-independent phenomenolog- 
ical scattering rate i/r. Note that the second term on 
the right-hand side of ( |95| ) is absent in Ref. |3(| because 
there the diamagnetic susceptibility of the electron gas, 
Xii c (q, 0), was implicitly taken to be zero. 

Finally, explicit expressions for the memory functions 
Af L (q, u>) and M T (q, u>) are obtained in the following 
way from Eq. (|70|): first, we approximately write 



M Jj (q,w) = (P Jj (q) 



1 



F L (q) 



(96) 



and similar for M T , i.e. we resort to the standard ap- 
proximation of replacing the projected by the full Liou- 
ville operator by setting Q sa 1 in the denominator. We 
thus assume that the huge amount of degrees of free- 
dom in the thermal bath and their extremely complex 
time evolution are completely dominating over the small 
subspace of observables Ao,Ai and its relatively well- 
controlled time evolution. 

In turn, the fluctuating longitudinal forces have the 
form 

V™ ^ q- 
and the transverse forces are 

F T (q) = -^E qx(q 2 xk) ^(q-k)p(k). (98) 

Here, U (q) is a random scattering potential. In the weak- 
disorder limit, we can perform the following decoupling, 
up to within corrections of higher than second order in 
the disorder potential: 
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[/(q-k)p(k) 



C 



[/(q-k'Wk')U W 



/0 (k'))(C/(q-k)C7(q-k')) . (99) 

disorder 



Since the system is homogeneous, we have k = k', and we arrive at the following expression for the longitudinal 
memory function: 



M^w) = ij>(q-k)) 2 [q ' (q 2 k)]2 C&(k jW ) x &(k, 



0) 



(100) 



One thus needs to calculate M L (q, lo) and C^ (q, to) via self-consistent solution of Eqs. J87| ) and ( |l00| ). The so 
determined Cr^q, ui) then serves as input for the transverse memory function 



M T ( q ,Lo) = l]T([/(q-k)) 2 [qX (q 4 X k)]2 C&Ck.aOxfcCk.O) , 



(101) 



which was obtained using the same decoupling approximation that led to Eq. (tlOOh for M (q,w). 



IV. MEMORY FUNCTION FORMALISM AND TCDFT FOR INHOMOGENEOUS SYSTEMS 

A. Formalism 



We now generalize the memory function formalism to systems that are inhomogeneous in one spatial direction, but 
still homogeneous in the plane perpendicular to it. The example we have in mind are quantum wells whose direction 
of growth is the z-axis. One can then in general no longer decouple longitudinal and transverse components of the 
correlation functions. A special case where this is still possible will be discussed in some detail later on. 

The generalization of Eq. (In]) for this inhomogeneous situation is 



k 



dz 1 ' 



u 5{z" - z) S ik - Q ik (q_\\,z,z") + Mik(q.\\,z,z",u>) C k j(q.\\,z",z',u) = S(z - z'), 



(102) 



where q|| is the in-plane wave vector. We again consider density fluctuations and current density as the only variables. 
Just like in the homogeneous case, by virtue of the continuity equation, all elements of the 2x2 memory function 
matrix Mi k are zero except M\\. As before, we can separate intrinsic and extrinsic contributions to the memory 
function in the weak-disorder limit, M\\ = M- m + M ox , and we combine the intrinsic part with £li k , defining 



^.(q||,z,z',w) = Q, ik (qn,z,z') + M in (q||, z, z',u))8a5 k x 



(103) 



Eq. (102) thus becomes 



l 

E 

fc=0 



dz' 1 



ujS(z - z) + M cx (q||,z,z ,uj) 5 ik - O, $ (q| | , z, z ,u) >C k j(q.\\, z ,z ,u) 



SijS(z-z') + 5 i0 dz" M cx (q\\,z, z" ,w)Cij(q\\, z" , z' ,uj) , i,j = 0, 1. 



(104) 



Similar to Sec. [II C, we will solve this set of equations by introducing suitable reference functions. However, the 
inhomogeneity of the system prevents us from using the same trick as for the homogeneous case, where we directly 
expressed the correlation functions in the presence of disorder in terms of the "clean" correlation functions, with their 
frequency argument lo replaced by u) + M cx . Now, by contrast, the memory function M ex (qi| , z, z', to) is no longer 
simply a number, but acts in conjunction with an integral operator, see Eq. (104). To deal with this difficulty, we 
first define a set of intermediate reference functions C^(q^,z, z' 7 uj) that satisfy the following coupled equations: 



fe=0 



I dz" \ cj S(z" -z) + M cx (q, | , z, z", u>)\ S ik - n&(q„ , z, z", lo) \C% (q, | , z", z' , lo) 

= -Sij S(z - z') , i, j = 0, 1 . 



(105) 
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In terms of these reference functions, the full correlation functions are given, combining Eqs. (104) and (105), through 
the following Dyson-type integral equation: 

dz\ Jdz 2 C%(q.\\,z,zi,u)M ex {q\\,zi,Z2,u)CQj(ci\\,Z2,z',Lj) . (106) 
We now use a similar trick to obtain the reference functions C R . Defining the "clean" response function via 



fdz" L 6{z" 



z)8 lk - ^(q||,z,z",c) \C c kj (q ll ,z",z',co) = -S tJ S(z - z') 



(107) 



(108) 



we get from Eqs. ( fL06| ) and ( |T07| ) : 

Cij(q\\,z,z',uj) = C?j(q.u,z,z',u)) 

+ X! dz i dz 2 Cifc (q||, 2, ^1,^) M cx (q||,zi,z 2 ,w) C^(q||,z 2 ,z',a;) . 

The density-density reference correlation function is explicitly given by 
Ca {<l\\,z,z',u)) = C% (q\\,z,z',uj) 

+ J dzi J dz2 CooCqih^^i^) M ex (qn,zi,z 2 ,uj) C^(q||, z 2 , z', u) 

+ j "dzi j dz 2 Coi(q||, z, z x , w) M ex (q||, z t , z 2 , u) C^ Q (q\\, z 2 , z' ,lu) . (109) 

At this point, it is convenient to introduce the following set of auxiliary functions, which will allow us later to write 
the memory function M ex in a more compact form (see below): 

(110) 
(111) 
(112) 

Again it is assumed t hat all extrinsic damping effects can be neglected for the static response function. In terms of 
these functions, Eq. (106) becomes for i = j = 0: 

®(q\\,z,z',U)) = <S> R {q\\,z,z',uj) 



$(q||,z,z',w) = J dz" C q((1\\,z,z",uj)x c ((1\\,z",z',0) 
$ fl ( q|| ,z,z',aO = jdz"C R {q^Z,z",Lj) X C {<l\\,z",z',Q) 
$ c (q||,z,z', W ) ee /dz"C c (q||,z,z",o;)x c (q||,z",z',0) 



dzi I dz 2 dz 3 $ R (q\\,z,Zi,<jj) X c (Q|h z i, z 2, 0) M ox (qn , z 2 , z 3 , u) $(qii, z 3 , z', u) . 



(113) 



In the example to be discussed in section [v|, only the case q|| = will be of interest. Since i n this case the continuity 
equation can be used to explicitly eliminate the current density in favor of the density, Eq. ( 109 ) can be written as 

$ fl -(0,z,z',w) = $ c (0,z,z',w) 

+ Jdz 1 Jdz 2 Jdz 3 <S> c (0,z,z 1 ,uj) x c (0,zi,z 2 ,0) M cx (0,z 2 ,z 3 ,uj) <S> R (0, z 3> z' ,uj) 

dz x X c (0 7 z 7 zi^) M cx(0,zi,g 2 ,^) [ w ^ ( o,5 2 , z »_ x (o,z 2 ,z',0)] . (114) 

Wn(zi)n(z 2 ) J- 



+ / dz\ \ dz 2 



The desired density-density response functions are then finally obtained using 

$(<l\\,z,z',u) = [x(q\\,z,z',w) - x c (q\\,z,z',0)]/uj 
$ fl (q||,z,z',w) = [x R (q\\,Z,z',cj) - x c (q\\,z,z',0)]/u 
$ c (q||,z,z',w) = [x c (q\\,z lZ ',uj)-x c (<l\\, z , z ',0)}/u 



(115) 
(116) 
(117) 
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To summarize: Eqs. (|109| )- (|117| ) allow one to express the interacting density-density response function of the system 
in the presence of intrinsic and extrinsic dissipation, x(Q|| > z i z \ w )i m terms of the interacting response function for 
the "clean" system, x c (q^, z, z',uj), i.e. including intrinsic dissipation alone. \ c is calculated, in principle exactly, 
using the framewor k of TCD FT o utlined in Sec. ||. Although admittedly somewhat frightening in appearance, the 
integral equations (113) and (114) involve only one-dimensional integrals, and their numerical solution is therefore 
quite manageable, as will be shown below. 



B. Memory functions for impurity and interface roughness scattering 



In general, a vector field V(r) is decomposed into longitudinal and transverse components as follows: V L (r) = 
-(1/4tt)V / dV[V' ■ V(r')]/|r- r'| and V T (r) = (1/4tt)V x V x J dVV(r')/|r - r'|. In our case, i.e. working in a 
mixed (q|| , z)-representation, it is convenient to define the operator D(qii,z) as 



D(q,|,z) 



qi 

«V 2 



(118) 



The longitudinal fluctuating forces are then given by 

1 



{ D (qih z 

and the transverse fluctuating forces are 

1 



D(qii,z) dz' 



P\\ e" 



-qu\z-z'\ 



(2tt) 2 2?,| 
P(P||^') D (P|| -%,*f)U{q\\ -V\\,z') } 



(119) 



F J (qn^) = - 



D(q||, z) X D(q||, z) X / dz' 



P\\ e " 



(2tt) 2 2<zh 

/S(pii^') d (ph -qii^'Mqi ~Pih z ' 



(120) 



One finds from Eqs. (11£) and (120) that in the homogeneous limit (no z dependence of the density fluctuations p 
and the scattering potential [/), the longitudinal and transverse forces only have in- plane components, given by the 
2-D versions of Eqs. ( |97| ) and (pq) . In the following, we shall limit the discussion of the inhomogeneous situation to 
a case of special interest, namely q|| =0. In that case, only the z-component of F L survives, and is given by 



PI I 



VW)J ( 27r ) 

The transverse force, in the same limit, acts in the x — y plane only: 



p(pn,z)V z U(-pn,z) 



d 2 p\\ 



Pll P{V\\,z)U{-p\\,z 



(121) 



(122) 



We n ote that, by symmetry, in the limit q|| = there is a n atural decoupling of the formalism outlined above (Sec. 
IV A) into separate sets of equations of the type ( 113)— (117 ) determining longitudinal and transverse response func- 
tions, respectively. In other words, L — T cross correlations are absent since the associated fluctuating forces are 
perpendicular to each other. The longitudinal and transverse memory functions are obtained in a quite straightfor- 



ward manner from Eqs. ( 
homogeneous case in Sec. 



119|) and (120), using the same approximate decoupling procedure that was used for the 



UK 



and 



M^(z,z',ij) - 



The result is: 

d 2 P\\ $>{p\\,z,z',uj) 



(2tt) 2 y/ n (z)n(z') 



V z V^(t/(- P ||,z)t/(-p||,0 



d 2 pu $(pi|, z, z',w) 



(2tt)2 ^/ n (z)n(z>) 

r 



=^ ([/(-pii,zM-pi 



(123) 



(124) 



The next step consists in finding explicit forms for the disorder- averaged random scattering potential U(p\\,z), 
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associated with some extrinsic damping mechanism. In 
the following, we shall focus on two examples specific to 
quantum wells: damping by charged impurities and by 
interface roughness. 

The potential associated with a single, statically scree- 
ned, positively charged impurity at position z\ is 



2vr 



-P|||z-Zl| 



(125) 



-On) v\\ 

where s(p\i) is the 2D dielectric function.0 The longitu- 
dinal memory function for charged-impurity scattering is 
thus 

&V\\ ®(P\\,z,z',u) 



Mj(z,z',uj) = 



''(p\\) y/n(z)n(z') 



=- J dz rii(z) 



xsign(z — z)sign(z' 



z)e 



-P\\\z- 



-pu\z'-z\ 



,(126) 



where rii(z) is the number of impurities per volume. 

Likewise, the longitudinal memory function associated 
with interface roughness is 



(2tt) 2 y/n(z)n(z') 



x\7 z X7 z , [S(z - zi)S{z' - zi) + S(z - z r )S{z' - z r )](,127) 

where U(p\\) is the random roughness scattering poten- 
tial, assumed for simplicity to be the same at the left 
and right interfaces, zi and z r . It is common to assume 
a Gaussian form for the autocorrelation function of the 
random interface roughness ,c3c3 which leads to 



<C/(P||) 2 ) =VAVe 



2 ^ 2 -P|i 7 ? / 4 



(128) 



Here, \i is the height of the potential step at the interface, 
and the correlation length rj and average roughness height 
A are controlled by material and growth conditions. In 
the presence of both impurity and roughness scattering, 
the memory functions M\ and Mr are additive (i.e., dif- 
ferent extrinsic scattering mechanisms are assumed to be 
uncorrelated) . 

Some practical complications arise from the fact that 
the memory functions explicitly depend on 3>(qii, z, z', u>) 
at all q||, not just q|| = 0. This means that $ should be 
calculated self-consistently from Eqs. (113) for all q||, 
which is a very demanding computational task. There- 
fore, as a first approximation, we ignore self-consistency 
and instead use the noninteracting $q( ( 1||j z, z' , ui) in 



(126) a nd ( 127), defined by replacing x and Xc with x KS 
in Eq. (115). The wave vector dependence of <&o is thus 
known analytically (see the explicit expression for % KS be- 
low) . This is expected to be a reasonable approximation 
as long as plasmon damping is not too strong. 



V. LINEWIDTH OF INTERSUBBAND 
PLASMONS IN A QUANTUM WELL 

In semiconductor quantum wells, the conduction band 
splits up into several subbands, and electrons (supplied 



e.g. by remote doping) can perform collective transitions 
between them. These so-called intersubband (ISB) plas- 
mons arp-currently of great experimental and theoretical 
interest J£3 being the basis of a variety of new devices op- 
erating in the terahertz,*egime, such as detectorsc3 and 
quantum cascade lasers .ca In designing these devices, the 
emphasis usually lies in covering a particular frequency 
range. However, often it is desirable that the transitions 
also have a narrow linewidth, to achieve better frequency 
resolution and larger peak absorption in detectors, and 
higher gain in lasers. The linewidth arises from a com- 
plicated interplay of a variety of scattering mechanisms, 
intrinsic (electron-electron and electron-phonon) as well 
as extrinsic ones (impurity, alloy-disorder and interface 
roughness) . Many aspects of this interplay are still not 
well understood, in particular the relative importance of 
the individual mechanisms .E2I 

To disentangle the various contributions to the ISB 
linewidth, it is helpful to consider a situation where 
some of them axe not effective. In a recent experiment, 
Williams et alJEl studied collective ISB transitions in an 
n-type 40 nm wide single GaAsZAlo.3Gao.7As quantum 
well, with Si doping centers 100 nm away from the well. 
Sharp transitions were found well below the LO phonon 
frequency of GaAs (35.6 meV), at a temperature of 2.3 
K. Thus, neither remote impurity nor phonon scattering 
are playing any significant role (nor is alloy-disorder scat- 
tering, as shown in Ref. [5l]). The linewidth is therefore 
expected to be dominated by bulk impurity and interface 
roughness scattering, while electronic many-body effects 
have traditionally been neglected. However, for high- 
quality samples such as the one used in the experiment 
discussed here, this iSpiio longer justified. 

In the experiment, till two parameters were controlled 
independently: the electronic sheet density N s (from 0.05 
to 1.3 x 10 11 cm -2 ), and a static electric field E perpendi- 
cular to the well which pushes the electrons against one 
of its edges. This provides an ideal tool to distinguish 
interface roughness from other damping effects. 

We describe ISB plasmons within a one-band-effective- 
mass approximation with parabolic subbands ,E3 which is 
a widely used and, for our purposes, sufficiently accurate 
method for GaAs/Al .3Gao.7As quantum wells. The case 
of the clean quantum well was treated in detail in Ref. 
p3| . The non-interacting response function reads 

/Voce OO 

XKs(q||,z,z',w) = 'J2'J2F flv (q\\,uj)(p^(z)(p fl (z') 

xtp„(z)p„(z') , (129) 



where 



/(e^ + fc^/2) 



q||k|| +a llll {q\\) + oj + ir\ 



q||k|| + a llv (q\\) - u> - ir) 



, (130) 
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FIG. 1: Dispersions of ISB (circles) and 2D plasmon 
(squares) frequencies f2(?ii) in a clean quantum well. UJ12 
is the difference between the two lowest bare subband levels. 
Inside the shaded regions, the plasmons are subject to strong 
Landau damping and rapidly rdie off (open symbols) . The 
experiment by Williams et 0,1x3 measures the ISB plasmon 
frequency and linewidth at = 0. 

®ij,u(q\\) = + ( v ~ / is the Fermi function at 

T = 0, and rj is a positive infinitesimal. e M and f^{z) are 
the Kohn-Sham energies and wave function (in LDA) of 
the quantum well. For the experimental range of N s , the 
system under study has 9 .bound levels, only the lowest 
being occupied (N occ — l).E2l 

We consider perturbations of the form v cx t.i(z,Lu) = 
Eq z, corresponding to monochromatic plane electro- 
magnetic waves of amplitude Eq polarized along the 
z-axis, the direction of growth of the quantum well. 
Having solved the response equation (|f0|), the pho- 
toabsorption cross section is then obtained as <j{lo) = 
— (8TTUj/Eoc)lmJdzzni(z,uj) and can be directly com- 
pared with data from photoabsorption measurements. 
(j(w) has a peak at the plasmon frequency f2 with 
linewidth (HWHM) T. 

In Fig. [l]we plot the dispersions f2(<7||) of the ISB and 
the intrasubband (or 2D) plasmon in the clean quantum 
well (N a = 1.0 x 10 11 cm" 2 ), calculated within ALDA, 
see Eq. (|l5|). The imaginary part of the Kohn-Sham re- 
sponse function x KS determines the regime of damping 
by single-particle excitations (Landau damping) , as indi- 
cated by the shaded region in Fig. [I]. Outside that region, 
in particular at small q> | , the plasmons are undamped in 
ALDA, for which / xc is frequency-independent and real. 
In reality, however, the absence of momentum conserva- 
tion in the z-direction and coupling via Coulomb inter- 
action opens the possibility of plasmon decay into more 
complicated excitations, such as multiple electron-hole 
pairs, even at g| — 0. To take this effect into account, 
one has to go beyond the ALDA and include dynamical 
xc effects. 

In Fig. U we show the electric field dependence of the 
ISB plasmon frequencies fl{q\\ — 0) for different values of 
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FIG. 2: ISB plasmon frequency Q, at q\\ = 0, versus electric 
field E. Top: experimental data from Ref. [23]. Middle: cal- 
culated results for the clean quantum well, using TDFT [Eq. 
(p"o[)]. Bottom: bare LDA subband splitting. The individual 
curves are associated with different electronic sheet densities 
(N s = 0.05, 0.1, 0.2, 0.3, 0.5, 0.7, 1.0, and 1.3 x 10 11 cm" 2 ). 
The lower N 3 , the steeper fl(E) around E — 0. 



N s . In the experimental data, built-in electric fields are 
subtracted, so that fl(E) exhibits a minimum for E = 
and rises quadratically for small fields. fi(E) increases 
most rapidly for the smallest N s , since higher electronic 
densities tend to screen the external electric field more 
efficiently. At the same time, the depolarization shift 
increases with N s . As a consequence, the curves of Q(E) 
for different _/V s are crossing each other. These features 
are very well reproduced by theory. Ignoring dynamical 
xc effects (i.e. using RPA) induces a 10% blueshift of 
f2, which then compares less favorably with experiment. 
To demonstrate the importance of including many-body 
effects in the response equation, we plot the bare subband 
spacings in the bottom panel of Fig. ||. The results are 
clearly qualitatively wrong: there is no crossing of fl(E) 
for different N s , in contradiction to experiment. 

Fig. |^ shows the ISB plasmon linewidth T(E), for dif- 
ferent N s . For small E, the experimental data again 
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FIG. 3: ISB plasmon linewidth F, at qn — 0, versus elec- 
tric field E. Top: experimental data from Ref. [23]. Middle: 
calculated results including extrinsic (impurity and interface 
roughness) and intrinsic (electron-electron interaction) damp- 
ing. Bottom: Y{E) for a clean quantum well (intrinsic damp- 
ing only). The individual curves correspond to different elec- 
tronic sheet densities (see Fig. ^). At E = 0, the lowest N a 
causes the smallest F. 



exhibit a quadratic behavior, and T(E) rises faster for 
smaller N s . For large negative E, T saturates around 
0.7 meV. For positive E (i.e., pointing in the direction 
of sample growth), T rises somewhat higher. The asym- 
metry of T(E) is likely to be due to slightly different 
roughnesses of the interfaces. 

The calculated T for the clean quantum well (intrin- 
sic damping through electron-electron interaction only), 
shown in the bottom panel of Fig. [| lies clearly be- 
low the experimental values, which is hardly surprising. 
However, it can be seen that these purely electronic ef- 
fects are far from negligible, at least for N s not too small, 
and provide an intrinsic lower limit to the linewidth of 
order 0.1 to 0.2 meV for N s ~ 10 11 cm~ 2 . 

The middle part of Fig. || shows T(E) calculated 
including electronic, impurity, and interface roughness 
damping, using the combined TCDFT and memory func- 



tion formalisms outlined above, see Eqs. (113-117). The 
results are now in very good agreement with experiment 
even away from flat band, as long as E is not too large. 
We now briefly discuss the details of how the various con- 
tributions to extrinsic damping were modelled. 

The presence of bulk impurities in the quantum well 
is mainly caused by segregation of donors from the lower 
delta-doped layer and diffusion along z during growth. 
We use the functional form n (z) — 1.33 



s/30n 



10 cm , proposed in Ref. 31 to explain in-plane mo- 



bility data, for the bulk impurity concentration in Eq. 
( |126| ). We also include scattering from the upper delta- 
doped layer (remote impurity density 4.8 x 10 11 cm -2 ). 

In Ref. |3l] the in-plane mobility was found to be dom- 
inated by bulk impurity scattering. By contrast, it turns 
out that neither bulk nor remote impurities contribute 
much to the linewidth. The behavior of T is instead dom- 
inated by interface roughness scattering and can in fact 
be qualitatively explained by it alone: Via (127), T de- 



pends on the product of density fluctuations at the edges, 
which, for E = 0, have largest amplitude for highest iV s . 
For finite E, electrons get pushed towards one edge, but 
less so for higher densities due to screening of the exter- 
nal field. T(E) thus rises more steeply for smaller N s , 
and the curves cross. 

W e take a roughness scattering potential of the form 
(128). The height of the potential step for our quan- 
tum well is fj, = 257.6 meV. The roughness parameters 
are chosen as n — 64. 4A and A = 4A, to give the best 
fit to experiment for the largest N s . Both 77 and A 
are in the characteristic range found by lattice imaging 
techniques .E3 

We also find that including electron-electron scattering 
does lead to a significant quantitative improvement for T, 
in particular for small E. 

For \E\ >j 1 mV/nm, the experimental linewidth sat- 
urates. This saturation can be understood as a nega- 
tive feedback effect, rela ted to the self-consistency of the 
memory functions (126) and (127). Roughly speaking, 



the plasmon linewidth comes from the imaginary part of 
M L , which in turn depends on the imaginary part of $. 
Broadening of the plasmon resonance means that Im<J> is 
peaked around Q. over some frequency range of width T. 
But, due to the constraint of the /-sum rule, increasing T 
means that the height of the peak of Im$ must decrease. 
This, in turn, limits the growth the memory function and 
rapidly saturates T. Neglect of the self-consistency of the 
memory function, as in our calculations, thus means that 
this saturation effect cannot be fully captured, as can be 
seen from Fig. |[ 



VI. CONCLUSION 

In this paper, we have dealt with a long-standing prob- 
lem in the many-body theory of extended systems: the 
calculation of collective electronic excitations and their 
associated linewidths in systems that are both inhomo- 
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geneous and weakly disordered, in the sense that the ran- 
dom potential can be treated as a perturbation. In any 
real extended system, collective excitations are subject 
to dissipation, causing the associated coherent, plasmon- 
like motion to decay into many individual, incoherent 
degrees of freedom associated with single-particle excita- 
tions. It has been common practice in the literature to 
describe these processes with phenomenological assump- 
tions of varying degrees of refinement. This paper, by 
contrast, presents a new formalism that allows one to 
calculate excitation energies and their lifetimes entirely 
from first principles. 

Our approach begins with the basic notion that there 
are two classes of mechanisms that are responsible for 
dissipation of collective electronic dynamics. The first 
class is intrinsic scattering, which occurs even in a "per- 
fect" material or device. Here we have in mind primarily 
systems where phonon scattering is inactive, so the only 
source of intrinsic dissipation are electronic many-body 
effects such as multiple particle-hole excitations. The sec- 
ond class of dissipation mechanisms is extrinsic in nature, 
such as scattering off impurities and disorder. 

Our treatment of intrinsic scattering relies on TDFT 
for the linear response. Fundamental existence theorems 
guarantee that TDFT describes the linear dynamics of 
interacting many-electron systems in principle exactly, 
including dissipation of collective degrees of freedom. In 
practice, however, the success of a TDFT approach relies 
on the approximations used for the linearized xc poten- 
tial. The most widely used approximation, the ALDA, 
has proved to be useful for calculating accurate excitation 
energies, but it produces linewidths that are strictly zero. 
Thus, a non-adiabatic description, which includes retar- 
dation, is required. A non-adiabatic dynamical density- 
functional approach which is local in space but nonlocal 
in time has to be formulated replacing the density with 
the current as basic variable (TCDFT). As a result, the 
linearized xc potential in the TCDFT response equation 
in general acquires a frequency dependence and an imag- 
inary part, leading to finite linewidths. 

To deal with extrinsic scattering, on the other hand, 
we make use of a powerful formal technique, the so- 
called memory function formalism. This approach can 
be traced back to the relaxation time approximation, but 



it replaces the simple phenomenological relaxation time 
t with the memory function M(q, uS), which is defined 
microscopically as a correlation function between fluctu- 
ating random forces. The memory function formalism 
is developed in the language of Kubo relaxation func- 
tions, which are, however, intimately connected to the 
(current)density response functions. 

The final step then consists in uniting the mem- 
ory function formalism with linear response theory in 
TCDFT. We thus arrive at a new, self-consistent theory 
that expresses the response function of an interacting sys- 
tem in the presence of both intrinsic and extrinsic damp- 
ing in terms of the "clean" interacting response function 
(which contains only intrinsic damping) and the memory 
function (which accounts only for extrinsic damping). 

We finally applied the theory to describe ISB plas- 
mons in a wide GaAs/Alo.3Gao.7As quantum well. Using 
reasonable values for the roughness parameters, we ob- 
tained quantitative agreement with the experimentally 
measured linewidth. But we also found that purely elec- 
tronic damping due to dynamical exchange and correla- 
tions makes non-negligible contributions to the linewidth, 
especially for high electronic densities, where the effect 
can be as high as a few tens of percents. 

A further remarkable outcome of this study is the phys- 
ical insight that the ISB plasmon linewidth is primar- 
ily controlled by interfacial roughness, and only weakly 
affected by the concentration of bulk impurities. The 
opposite is true for the in-plane mobility, which is pri- 
marily controlled by bulk impurities.Ea Thus, the corre- 
lation between ISB plasmon linewidth and in-plane mo- 
bility is rather weak, which is physically understandable 
since currents are flowing perpendicular to the quantum 
well in the former case, and parallel to it in the latter. 
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